Self-similar disk packings as model spatial scale-free networks 
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The network of contacts in space-filling disk packings, such as the Apollonian packing, are ex- 
amined. These networks provide an interesting example of spatial scale-free networks, where the 
topology reflects the broad distribution of disk areas. A wide variety of topological and spatial 
properties of these systems are characterized. Their potential as models for networks of connected 
minima on energy landscapes is discussed. 
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I. INTRODUCTION 

Since the seminal paper of Watts and Strogatz on 
small- world networks [1], there has been a surge of in- 
terest in complex networks, both to characterize real- 
world networks and to generate network models to de- 
scribe their properties [2-7]. The systems analysed in 
this way have spanned an impressive range of fields, in- 
cluding astrophysics [8] , geophysics [9] , information tech- 
nology [10], biochemistry [11, 12], ecology [13] and sociol- 
ogy [14] . Initially, the focus was on relatively basic topo- 
logical properties of these networks, such as the average 
separation between nodes and the clustering coefficient to 
test whether they behaved like the Watts-Strogatz small- 
world networks [1], or the degree distribution to see if 
they could be classified as scale- free networks [15]. 

As the field has progressed, however, the emphasis has 
shifted away from these basic classifications to increas- 
ingly detailed characterization of the networks. For ex- 
ample, on a topological level, there has been much re- 
cent interest in both the correlations [16] and community 
structure [17] within a network. 

There has also been increasing interest in how the 
medium in which a network is embedded influences the 
network properties. For spatial networks this can of- 
ten lead to some kind of geographical localization [18]. 
For example, in social networks, acquaintances are more 
likely to share the same neighbourhood, and for the in- 
ternet there is obviously a greater cost associated with 
making longer physical connections [19]. To model these 
kinds of effects there have been a number of studies in 
which the preferential attachment rule that leads to scale- 
free networks [15] has been altered to include an addi- 
tional distance dependence in the attachment probabil- 
ity [20-23]. Typically, this leads to some crossover away 
from scale-free behaviour when the distance constraint is 
sufficiently strong. 

A different approach to understanding the interplay 
of geography and topology has been to consider ways in 
which a scale- free network can be embedded in Euclidean 
space [24-27]. In most of these spatial scale- free net- 
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works, the nodes are distributed homogeneously in space 
[24-26] . The heterogeneity that leads to the scale-free be- 
haviour instead comes from the node dependence of the 
range of interactions, i.e. high degree nodes have con- 
nections to nodes that lie within a larger neighbourhood 
of the node. The model of Herrmann et a/., however, 
shows the converse behaviour [27]. Each node has the 
same interaction range; instead the scale-free behaviour 
is driven by an inhomogeneous density distribution with 
high-degree nodes associated with regions of high node 
density. 

Our interest in spatial scale-free networks comes from 
recent work characterizing the connectivity of the config- 
uration space of atomic clusters [28, 29]. Configuration 
space can be divided up into basins of attraction sur- 
rounding each of the minima on the potential energy sur- 
face of the clusters [30]. This then allows a network de- 
scription of the potential energy surface where the nodes 
correspond to the minima, and two minima are linked 
if there is a transition state valley directly connecting 
them. All links are therefore between adjacent basins of 
attraction. Intriguingly, this "energy landscape" network 
was found to be scale free. Since that initial study, the 
configuration space of some polypeptide chains has also 
been found to have a scale- free connectivity [31]. 

This scale-free behaviour cannot be explained by 
the usual preferential attachment approach [15] because 
these networks are static, and are just determined by the 
potential for the system. Neither are the spatial scale- 
free models described above much help, because they are 
not contact networks between spatially adjacent regions. 
For example, if one were to associate each point in Eu- 
clidean space of these models with the nearest node the 
network of contacts between the resulting cells would not 
be scale-free. Instead, the scale-free behaviour of these 
spatial networks arises precisely because there are more 
long-range connections between non-adjacent nodes. 

To try to further understand the energy landscape net- 
works a different approach must be taken. In Ref. [28] it 
was suggested that the scale-free behaviour might reflect 
differences in the basin areas, with the deeper minima 
having large basins of attraction [32] with many connec- 
tions to the smaller basins surrounding them. For this 
to lead to a scale-free topology, one would imagine that 
the basins have to be arranged in some kind of hierarchi- 
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FIG. 1: An Apollonian packing of disks within a circle. 



cal fashion with basins at each level being surrounded by 
successively smaller basins. 

Space-filling disk packings, such as the Apollonian 
packing depicted in Figure 1, have just such features. 
In this paper, we examine the contact networks for such 
packings to determine whether they might provide a use- 
ful model for the energy landscape networks. In the final 
stages of the preparation of this work, Andrade et al. in- 
dependently introduced the idea of Apollonian networks 
[33]. In that work, only a brief characterization of the 
topology of the two-dimensional (2D) Apollonian packing 
was given, before the emphasis switched to the behaviour 
of dynamical processes on these networks. Here, we pro- 
vide a much more detailed characterization of the topol- 
ogy of the 2D network (Section II A), and also analyse 
the networks associated with other self-similar circle and 
hypersphere packings (Section II B). Furthermore, as our 
aim is to provide a model to help understand the energy 
landscape networks, a particular emphasis is the relation- 
ship between the topological properties of the networks 
and the spatial properties of the packings (Section III). 

II. TOPOLOGICAL PROPERTIES 

A. 2D Apollonian networks 

To produce an Apollonian packing, we start with an 
initial array of touching disks, the interstices of which 
are curvilinear triangles. In the first generation disks are 
added inside each interstice in the initial configuration, 
such that these disks touch each of the disks bounding 
the curvilinear triangles. The positions and radii of these 
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FIG. 2: (colour online) The development of the 2D Apollonian 
network inside the interstice between three mutually touching 
disks, as the number of generations increases. In each picture, 
the network is overlaid on the underlying packing. 

disks can easily be calculated using either the Soddy for- 
mula [34] or by applying circular inversions. Of course, 
these added disks cannot fill all of the space in the inter- 
stices, but instead give rise to three smaller interstices. In 
the second generation, further disks are added inside all 
of these new interstices, which again touch the surround- 
ing disks. This process is then repeated for successive 
generations. If we denote the number of generations by 
t, where t = corresponds to the initial configuration, as 
t — > oo the space-filling Apollonian packing is obtained. 

There are two initial configurations that are commonly 
used. The first was used for Figure 1, and has three 
mutually-touching disks all inside and in contact with 
a larger circle. This configuration has the useful feature 
that all of the initial disks are inside a curvilinear triangle 
formed by the other three circles, so that their initial en- 
vironment is equivalent to that for any of the subsequent 
disks immediately after generation. This configuration is 
therefore more convenient when analytically deriving the 
properties of the packing, because for the most part the 
same formulae apply to the initial disks as to all other 
disks. However, the spatial nature of the connections to 
the bounding circle are more ill-defined 

The second common initial configuration is just to have 
three mutually touching disks as in Figure 2, and in sub- 
sequent generations to progressively fill the single curvi- 
linear triangle in the initial configuration. The proper- 
ties of the initial disks do not follow exactly the same 
trends as the subsequent disks, but it has the advantage 
that none of the disks touch the interior boundary of 
another circle. The numerical results presented will typi- 
cally be for this initial triangular configuration. However, 
we should emphasize that the networks associated with 
these two initial configurations have the same topological 
properties, except for the initial disks. 

Figure 2 illustrates how the Apollonian packing can be 
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used as a basis for a network, where each disk is a node in 
the network and nodes are connected if the corresponding 
disks are in contact. We shall call this contact network 
an Apollonian network. Of course, it has an infinite num- 
ber of nodes. However, we shall generally consider the 
network properties after a finite number of generations 
in the development of the complete Apollonian network. 
From our analytical results we can quickly obtain the 
properties of the complete network by taking the limit 
of large t. However, the numerical results are necessarily 
limited to finite-sized networks. 

Figure 2 also shows how the network evolves with the 
addition of new nodes at each generation. For each new 
disk added, three new interstices in the packing are cre- 
ated, that will be filled in the next generation. Equiva- 
lently, for each new node added, three new triangles are 
created in the network, into which nodes will be inserted 
in the next generation. Therefore, 

An v (t) = n v (t) -n v (t - 1) = 3An v (t - 1) t > 1, (1) 

where n v is the number of nodes in the network. 

For the Apollonian packing of a circle, n„(0) = 4, and 
An v (l) = 4, if we treat the bounding circle on the same 
footing as the other three initial disks. It follows that 

An„(i) = 43*" 1 and n v {t) = 2 (3* + l) (2) 

The addition of each new node leads to three new edges. 
Therefore, 



An e (t) = n e (t) - n e (t - 1) = 3An v (t) =4 3* 



(3) 



where n e is the number of edges in the network. As 
n e (0) = 6 



n e (t) = 2 3 t+1 



(4) 



p, the fraction of all the possible pairs of nodes that are 
actually connected, is given by 



2n e 



3 t+i 



n v (n v -l) (3* + 1)(3* + 1/2) 
« gt~r for large t. (5) 

Therefore, the Apollonian network becomes increasingly 
sparse as its size increases. By contrast, the average de- 
gree tends to a limiting value: 



(k) 



2n e _ 2 3* +1 
~ ~ 3* + 1 
6 for large t. 



(6) 



Surrounding each node i are fcj empty (i.e. not enclos- 
ing any nodes) triangles. As new nodes are added at the 
centre of all these triangles, ki new connections will be 
created. Therefore, at each step the degree of a node 
doubles, i.e. 



h{t + l) = 2ki(t). 



(7) 



Such a rule expresses a preferential attachment [15]. The 
number of new connections is linearly proportional to the 
degree. 

If i Cj i is the step at which a node i is created, k{t c ^) = 3 
and hence 

h(t) = 32*- t '-*. (8) 

Therefore, k can take a series of discrete values up to 
tax = 32*. The fraction of the other nodes that the 
node with maximum degree connects to is given by 



3 2* 



n v - 1 



23* + 1 
t-i 



for large t 



(9) 



and is a decreasing fraction of the total as the size of the 
network increases. It follows that the degree distribution 
is given by 



for k = 3 2*,i r 







p(k) 



{ MO) = 2 

n v (t) ~ 3* + 1 

^M = ^rl forfc = 32*-st c >i 

n v (t) 3* + l " 
otherwise 



and that the cumulative degree distribution is 

3*- + 1 



Pcum(k) 



3* + 1 ' 



(10) 



(11) 



Substituting for t c in this expression using t c = t — 
log(fc/3)/log2 gives 



Pcum(^) — 



3* (ft/3) - 



-log 3/ log 2 



3* + 1 

log 3 
log 2 



for large t. (12) 



For a continuous degree distribution p{k) ~ fc~ 7 , 
Pcum(^) ~ k 1 ^ 1 . Therefore, the Apollonian network is 
scale free and the exponent of the degree distribution is 



7=1 + 



log 3 
log 2 



2.585, 



(13) 



as already noted in Ref. 33. Apollonian networks hence 
provide a new model for spatial scale-free networks. Im- 
portantly, in contrast to other two-dimensional spatial 
scale- free networks [20-27], the Apollonian network can 
be embedded in a plane without any edges crossing. In 
Aste et aVs classification scheme, they hence have a 
genus of zero [35] . 

Another important property of the network is the clus- 
tering, which provides a measure of the local structure 
within the network. The clustering coefficient of node i 
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is the probability that a pair of neighbours of i are them- 
selves connected. 



Ci{t) = 



2n< 



h(t){h(t) - \y 



(14) 



where n^°™ is the number of connections between the 
neighbours of i. At each stage a ring of new connections 
passing through all the nodes connected to i is generated. 
Therefore, 



\t)= ]T fc i (t') = 3(2*- t -+ 1 -l), (15) 

t'=t Cii 



and 



Ci(t) 



2 (2*-*e.« + l _ 1) 

2*-*.,, (32*-*^ - 1) 
1 4 



32 t-tc,i-2 



for f > t C)i . (16) 



Therefore, the clustering coefficient of a vertex shows the 
same inverse proportionality to the degree as has been 
obtained previously for other deterministic scale-free net- 
works [36-38]. This feature has been taken to be a sig- 
nature of a hierarchical structure to the network [36, 39], 
but recently has been shown to partially reflect disas- 
sortative correlations [40]. In the current networks this 
feature can also be interpreted in terms of spatial local- 
ization. For a low-degree node its neighbourhood only 
occupies a small local region in the packing, and thus 
would be expected to have strong clustering. By contrast, 
high-degree nodes have a more global character and are 
connected to well-separated parts of the packing, and so 
have low clustering. 

The clustering coefficient for the whole graph can be 
defined in two ways. The first is a generalization of Eq. 
(14) to the whole graph, and is the probability that any 
pair of nodes with a common neighbour are themselves 
connected. Thus, 



2 EX 



(17) 



The second definition of C is as the average value of the 
local clustering coefficient, i.e. 



rt... t—^i 



(18) 



The difference between these two definitions is the rel- 
ative weight given to nodes with different degree. High 
degree nodes make a larger contribution to C\ because 
there will be more pairs of nodes that have a high-degree 
node as a common neighbour, whereas all nodes con- 
tribute equally to C 2 - Typically, C\ < C2 because, as 
is the case here, higher degree nodes tend to have lower 
values of &. 
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FIG. 3: The dependence of the vertex betweenness on degree 
for the Apollonian network with t = 8. There are points 
representing each disk and the line connects the average values 
for a given k. 



Substituting in and rearranging gives 

3 t+1 - 1 



C\ = 



12 (2 2 *" 1 - 3*" 1 ) 
t 



(19) 



The clustering coefficient goes down as the size of the 
Apollonian network increases. However when one com- 
pares C to that for an Erdos-Renyi random graph [41, 42] 
(Cer = p) , one obtains 



Ci 



ER 



ill 

6 \2 



(20) 



That is the Apollonian networks become increasingly 
more clustered than a random graph, as their size in- 
creases. C2 can be evaluated numerically. As shown in 
Ref. 33, C2 tends to a constant at large t of value 0.828. 

In Ref. 33, they also calculated the behaviour of i ave , 
the average number of steps on the shortest path between 
any two nodes. Z avo showed a small-world behaviour, scal- 
ing sub-logarithmically with network size [33]. The im- 
portant role played by the larger disks in mediating these 
short paths is illustrated in Fig. 3. The vertex between- 
ness of a node is defined as the fraction of all the shortest 
paths that pass through that node. For t = 8, 40% of 
these paths pass through the central disk in the pack- 
ing (Fig. 2). Furthermore, the dependence of the vertex 
betweenness on the degree is not far from a power-law. 
This type of behaviour is common for scale- free networks 

Correlations in networks, particularly with respect to 
the degree, have been the subject of increasing interest 
[16]. This is partly because the behaviour of models de- 
fined on such networks have been often found to depend 
sensitively not only on the degree distribution, but also 
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on how correlated the networks are [44-47]. For exam- 
ple, k nn (k), the average degree of the neighbours of nodes 
with degree k, should be independent of k for an uncor- 
rected network. 

We can calculate fc nn for the Apollonian network using 
Eq. (7) to work out how many connections are made at 
a particular step to nodes with a particular degree. Ex- 
cept for the initial disks, no disks created in the same 
generation, i.e. with the same degree, will be connected. 
All connections to nodes with higher degree are made at 
the generation step, and then connections to lower degree 
nodes are made at each subsequent step. This leads to 
the expression 



k n n{k) — 



1 



An v (t c )k(t c ,t) 



t c =t c -i 

An v (t' c )k(t' c ,t c -l)k(t' c ,t) + 

t' c =o 

t' c =t \ 

An v (t c )k(t c ,t' c -l)k{t' c ,t)) (21) 

t' c = *c + l / 

for k = 3 2*~* c and where k(t c , t) is the degree of a node 
at generation t that was created at generation t c . The 
first sum corresponds to the connections made to nodes 
with higher degree (i.e. t' c < t c ) when the node was cre- 
ated at t Cl and the second sum to the connections made 
to the current lowest degree node at each step t' c > t c . 
After substitution and evaluation of the sums, the above 
expression simplifies to 



knn{k) 



"5 



6 + 



(22) 



After the initial generation step k nn for a node increases 
linearly with age. 

Writing the above equation in terms of k gives 



M*0 = 9 ( - 



log(4/3) 
log 2 



„ 31og(fc/3) ,„„, 

6 + (23) 



k n n(k) is roughly a power law function of k with expo- 
nent —0.415. The exponent is negative implying that the 
network is disassortative, i.e. nodes are more likely to be 
linked to nodes with dissimilar degree. When normalized 
by the expected value of k for an uncorrelated network 



^uncorr 



(k) \z 



(24) 



k nn (k) has a universal form independent of network size 
for large t and small k. Namely, 



knn(k) 
iUuncorr 



log(4/3) 
3 fk\ l°g2 



2 V 3 



(25) 
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FIG. 4: Comparison of fc n n/fcn£ corr for the 2D Apollonian net- 
work to that for randomized versions of the networks. Lines 
corresponding to different values of t (t = 7-15 and 7-8 for 
the Apollonian, and randomized networks, respectively) have 
been plotted to emphasize the universal form of this function 
for the Apollonian network. 



This is illustrated in Figure 4. The upward curvature 
away from this power-law form at large k is caused by 
the logarithmic term in k in Equation (22). 

It has been shown that disassortativity can often arise 
in networks where self-connections and multiple edges 
are excluded [48]. Therefore, k nn was compared to that 
for random networks with the same degree distribu- 
tion, which were prepared using the switching algorithm 
[16, 49]. The randomized networks also show disassor- 
tativity, but to a somewhat lesser degree (Fig. 4). The 
additional disassortativity arises because connections to 
the nodes with the same degree cannot occur in the Apol- 
lonian network (except for the initial disks). 

In particular, as the nodes with k = 3 are only con- 
nected to higher degree nodes, fc nn (3) is significantly 
higher than that for the randomized graph (Figure 4). 
By contrast, fc nn for the rest of the network is lower than 
that for the randomized graphs. This is also because 
of the lack of same-degree connections, as this gives the 
higher degree nodes many more connections to the most 
numerous k = 3 nodes than for the randomized graphs. 

An assortativity coefficient r has been proposed that 
measures the degree of (dis)assortativity of a property 
[50], and is defined as 



r = 



{st) e - (s) e (t)e 
e,assort (s} e {t) t 



(26) 



where s and t correspond to the property of interest at ei- 
ther end of an edge, e denotes that the averages are over 
all edges and assort that the average is for a perfectly 
assortative network, r is therefore a measure of the cor- 
relations in the property compared to that for a perfectly 
assortative network. Disassortative networks have r < 0. 
It follows that the assortativity coefficient for the de- 
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FIG. 5: Assortativity coefficients, and ta, as a function 
of the number of generations t. The solid lines are for the 
2D Apollonian network, and the dashed lines for randomized 
versions of these networks with the same degree distribution. 



gree is given by 



2\2 



Tk = 



(fc )(fc 2 fc nn (fc))-(fc 2 ) 
(fc)(fc 3 ) - (fc 2 ) 2 ' 



(27) 



Expressions for the quantities in the above equation can 
be relatively easily obtained from the degree distribution. 
This gives 



Tk = 



3* (2t-3 + 4(|) t )-2 2 ' (2 -(f)' 



3*2* 
10 /2 



(e -(!)*) - 22t ( 2 -(!)*) 

for large t. (28) 



3 V3 



Tk is always negative, indicating disassortativity. How- 
ever, its magnitude goes to zero as the size of the network 
increases (Figure 5). This may seem surprising since k nn 
always has a negative slope and has an effective func- 
tional form that is independent of size (Eq. (25)). How- 
ever, the convergence of r to zero is simply because the 
denominator corresponding to the correlations in a per- 
fectly assortative network scales more rapidly with size 
than the numerator, i.e. as 6* compared to —4* (Eq. (28)). 

We also looked at the community structure of this net- 
work using the algorithm of Girvan and Newman [17]. It 
works by starting with the complete network and at each 
step removing the edge that has the maximum edge be- 
tweenness, where this quantity is recalculated after the 
removal of every edge and is defined as the fraction of all 
the shortest paths that pass through an edge. If there 
is more than one edge with the same maximum edge be- 
tweenness, they are all removed at the same step. Thus 
the network is progressively divided into n v communi- 
ties. To decide which division of the network represents 
the best choice, the modularity Q is calculated at each 




FIG. 6: The best division of the 2D Apollonian packing into 
communities for t = 6. The disks in the different communities 
have been displaced with respect to each other for clarity. The 
modularity Q = 0.5938. 



step, and the division of the system with maximum Q 
is considered to be the best [51]. Q is denned as the 
fraction of edges that are within the communities com- 
pared to that expected for a random graph with the same 
degree distribution. 

The best division of the packing into communities is 
shown in Figure 6 for t = 6 and has Q = 0.5938. This 
value is comparable to some of the higher values found 
for networks considered previously [51, 52]. Interest- 
ingly, the Apollonian network's combination of commu- 
nity structure and assortativity is in contrast to most of 
the other networks with high Q which tend to be also 
strongly assortative. As expected the communities are 
spatially localized. As the algorithm only used topolog- 
ical information, this result implies that the spatial em- 
bedding of the network is clearly reflected in its topology. 

The algorithm for detecting communities described 
above cannot actually break the threefold symmetry of 
the packing. The best threefold symmetric division has 
the central disk as its own community. However, by as- 
signing this disk to one of the adjacent communities, Q 
improves from 0.5872 to 0.5938. We have also applied 
the faster algorithm described in Ref. 52, however slightly 
lower values of Q were obtained. 

It is interesting to examine how the network is progres- 
sively broken into separate communities as more edges 
are removed. This can be represented by a dendrogram 
as in Figure 7, which shows the number of communities 
and their relationship at each step in the algorithm. Ev- 
ery time two sets of disks become disconnected, their cor- 
responding lines split. The first 53 sets of edges removed 
only break the packing into two communities; instead the 
effect is to make the contact matrix sparser. Further re- 
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FIG. 7: A dendrogram showing the progressive division of 
the 2D Apollonian packing into communities for t = 6 as 
those sets of edges with largest edge betweenness are removed. 
At right no edges have been removed and there is a single 
community of size n v , whereas at left all edges have been 
removed and there are n v communities of size 1 . The number 
of steps in which sets of edges have been removed increases 
linearly from right to left. In total 84 sets of edges were 
removed. The three-fold symmetry of the packing is retained 
throughout this process. The vertical dashed line indicates 
the set of communities with largest Q. The inset shows how 
Q varies as the network is broken up into communities. 



moval of sets of edges then relatively quickly breaks the 
packing into a large number of communities. 

It is evident that Q has quite a broad maximum as 
a function of the number of communities. Q is greater 
then 0.3 when there are between 6 and 88 communities. 
Given the self-similar nature of the packing, one would 
not expect there to be a strongly preferred size for the 
communities. Similarly in the deterministic scale-free 
networks the hierarchical modularity means that there 
is no clearly preferred size for the modules at which the 
modularity is significantly enhanced [39] . 



B. Other self-similar packings 

The 2D Apollonian network is only one example of a 
space-filling self-similar packing. In a similar way to the 
last section, a detailed characterization of the topology of 
contact networks associated with other self-similar pack- 
ings could be derived. However, here we do not wish to 
give such a comprehensive account, but to illustrate how 
some of the key features of these networks, particularly 
the exponent of the degree distribution, depend on the 
nature of the packing. 

Firstly we shall examine higher-dimensional Apollo- 
nian packings. The initial configuration that is directly 
equivalent to Figure 1 is to have touching hyperspheres 
at the corners of a d-dimensional simplex that is enclosed 
within and touching a larger hypersphere. The analysis 



of the last section is relatively easy to generalize to these 
cases. 

As, now n v (0) = An„(l) = (d + 2) and An v (t) = 
(d + l)An„(t - 1) for t > 1, it follows that 



i-i 



and 



An v (t) = (d + 2)(d+l) 



n v (t)= ( -l±^((d+l) t + d-l). 



(29) 



(30) 



The higher dimensional equivalent of Eq. (7) is not so 
useful for calculating the degree distribution, for exam- 
ple for a 3-dimensional Apollonian packing fc,(£ + 1) = 
3ki(t) — 4. Instead, an alternative approach has to be 
used. Each new neighbour of a node i creates d new 
e?-simplices involving i. In the next generation these d- 
simplices will be the sites for new nodes that are also 
neighbours of i. Therefore, 



Afci(i) = fe(t) -h(t- 1) 

= dAk t (t - 1) t > t Cti + 1. 

As ki{t c ^) = d + 1 and Ak t (t c ^ + 1) = d + 1, 

Aki(t) = (d + ljd*-**' -1 



and 



d+1 



- (d*-*^ + d - 2) . 



(31) 



(32) 



(33) 



By an equivalent analysis to that for d = 2 one can show 
that Pcum(k) follows a power-law for large t where 



7=1 + 



logjd + 1) 
logd 



(34) 



Hence, the Apollonian networks associated with higher- 
dimensional packings are also scale- free networks. The 
exponent 7 decreases as the dimension of the Apollonian 
packing increases, tending to two in the limit of large d. 
This is noteworthy since the value of 7 can have signifi- 
cant effects on network properties [53] . 

By physical arguments it is easy to see that these 
higher-dimensional Apollonian networks will have very 
similar topological properties to the two-dimensional case 
that we have studied in detail. The networks will again be 
disassortative with respect to degree because of the lack 
of connections between nodes with the same degree. The 
hierarchical structure and the more localized character of 
the connections involving low-degree nodes will lead to a 
strong dependence of the clustering coefficient on degree. 
This spatial localization will also lead to strong commu- 
nity structure. Larger hyperspheres are also more likely 
to have a larger degree. 

To illustrate how these conjectures can be backed up 
analytically here we derive a general expression for the 
local clustering coefficient. We first need to calculate the 
number of connections between the neighbours of a node. 
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1=2 



FIG. 8: A space-filling bearing between two half planes. The 
pattern can be repeated to the left and to the right ad in- 
finitum. In the numenclature of Refs. 54-56 this packing is 
from Family 1 and has base loop size 4 and n = m — 2. The 
half-planes at the top and bottom can be considered as disks 
with infinite radius. 



FIG. 9: (colour online) The development of the contact net- 
work for a subset of the packing in Fig. 8 initially consisting 
of four touching disks, one of which has infinite radius. The 
edges in the network connect the centres of the disks that 
touch. Those edges involving the infinitely large disk are par- 
allel, but to emphasise that they do actually connect to the 
same node, a box has been added that enclosed the ends of 
these edges. 



On generation a node is surrounded by a d-dimensional 
simplex, and then at subsequent steps every new neigh- 
bour of a node contributes d new connections to nf on . 
Hence, 



d(d+l) 



2 

d+ 1 
d-l 



Y dAki(t) 



t'=t c +i 



rfW1 + ^ + l)(d-3) 



2(d-l). 



Substituting into Eq. 14 and taking limits gives: 



Ci(t) 



2d 



for t ^> t c , 



(36) 



Again, the local clustering coefficient is inversely propor- 
tional to degree. 

In the Apollonian packing of disks, the smallest loops 
in the contact network have size three, i.e. they are trian- 
gles. However, space-filling packings of disks are possible, 
where the smallest loops in the contact network are poly- 
gons with more than three sides. Examples, where the 
loops all have an even number of sides are of particu- 
lar interest, since they can act as space-filling bearings, 
where all the disks can rotate at the same time without 
slip [55, 57]. An example of a space-filling bearings with 
'base loop size' 4 and l n — m\ is shown in Figure 8. The 
procedures to construct such packings are more complex 
than for Apollonian packings and are described in detail 
in Refs. 54-56 

Figure 9 illustrates how the contact network associated 
with this packing develops for a subset initially consisting 
of four touching disks. At each stage to + 1 new nodes 
are produced for each empty quadrilateral, dividing the 
quadrilateral into a further to + 2 new quadrilaterals to 
which new nodes will be added at the next generation. 

Assuming that there is one quadrilateral initially (i = 
0), the number of empty quadrilaterals after step t is 
(m + 2)*. Hence, 



An v {t) = (to + l)(m + 2) 



t-i 



(37) 



and 



n' v [t) = n v (t) - n v {0) = (to + 2)' - 1. 



(38) 



We exclude the initial boundary disks in the above be- 
cause they have slightly different properties from the rest 
of the disks in the iterative scheme. Besides, for suffi- 
ciently large t their contribution is negligible. 

Aside from the first step after a node is created, a node 
only gains new connections at every other step, because 
the new connections are added across alternating diago- 
nals of the quadrilaterals. 

ki(t + 2) = (m + 2)ki(t). (39) 

Therefore, 

ki(t) = 2(m + 2) {t ~ t ^ )/2 for t - t C}i even. (40) 
In the same way as before this leads to 

(to + 2)*(£:/2)- 2 - 1 



Pcum 



(to + 2)* - 1 



for t large. 



(41) 



Hence 7 = —3, independent of to. Unlike the Apollonian 
packings, An v and fcj increase by the same factor to + 
2, and so this cancels. The power of two in the above 
equation arises, because this increase in fcj occurs only at 
every other step. The situations is more complicated for 
the case m 7^ n, because there are two sub-populations of 
disks, and the degree distribution of each sub-population 
obeys its own power law. 

We should note that the contact networks for these 
space-filling bearings have no triangles, and so the clus- 
tering coefficient defined by Eq. (17) is zero. However, 
generalized clustering coefficients probing higher-order 
loops have been proposed. Clearly for these space-filling 
bearings the number of 'squares' (loops of length 4) will 
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be significantly higher than for a random network, al- 
though we have not sought to quantify this. 

The results in this section illustrate that the contact 
networks associated with other space-filling disk and hy- 
persphere packings are also scale-free, but that the expo- 
nent of the degree distribution is not a universal constant 
but depends on the nature of the packing. We could have 
also considered other examples, such as scale-free bear- 
ings with base loop size greater than four [57] and non- 
Apollonian packings of spheres [58, 59], and it is likely 
that these again show somewhat different behaviour. 



III. SPATIAL PROPERTIES 

The 2D Apollonian packing is a well-known example 
of a fractal [60] , and has many of the typical fractal prop- 
erties. For example, inside every curvilinear triangle no 
matter how small the same pattern of disk packing reoc- 
curs, i.e. it is self-similar. Similarly, the estimated total 
length of the circumferences of all the circles continues to 
increase as the resolution of the measurement increase. 
One of the most important properties of such a packing 
is its fractal dimension. To understand this quantity we 
need to define more carefully the set for which we wish 
to know the dimension. In the packing the disks are 
all considered to be open, that is the set of points associ- 
ated with a disk Di contains all the points inside the disk 
boundary, but not the boundary itself. The residual set 
R is then the points that are not part of any of the open 
disks in the packing, or more formally R = U — \J°1 1 Di 
where U is the set that is being packed, dp, the fractal 
dimension of i?, is the quantity of interest. 

dp must obey 1 < dp < 2. The upper bound is obvi- 
ous, because, by virtue of the space-filling nature of the 
Apollonian packing, R must have zero area. The lower 
bound follows from the fact that r « = 00 1 where 

Ti is the radius of disk i.e. the total length of the 
boundaries is infinite [61]. This result can most easily 
be visualized by projecting the boundaries of each disk 
onto a diameter of the circle bounding the region that is 
being packed. Points on this diameter are projected onto 
infinitely often. These bounds imply the dimension of R 
must be fractional, and hence the packing is fractal. 

So far, no analytic formula for the value of the frac- 
tal dimension for the 2D Apollonian packing has been 
obtained, but instead its numerical value has beens es- 
timated with increasing precision [54, 62-66]. Its value 
is 1.3057. It has been suggested that the fractal dimen- 
sion of the Apollonian packing is the minimum for any 
space-filling disk packing [67], because at each step of the 
generating process the disk with maximum possible ra- 
dius is inserted into each curvilinear triangle, thus max- 
imizing the area of the region that must be outside of 
the residual set. The fractal dimensions found for 2D 
space-filling bearings are consistent with this assertion 
[54, 57]; all have larger values than that for the Apollo- 
nian packing with the largest found being 1.803. In fact, 



as pointed out by Melzak, it is easy to generate a disk 
packing with dimension arbitrarily close to 2 [67] . If each 
disk in the Apollonian packing in Figure 1 is replaced by 
a suitably scaled image of the whole Apollonian packing, 
a new packing is obtained with higher fractal dimension. 
If this process is repeated ad infinitum, a disk packing 
with a fractal dimension of 2 is eventually obtained. 

For space-filling packings of d-dimcnsional hyper- 
spheres, there are similar limits for the fractal dimen- 
sion, namely d— 1 < dp < d. The only calculations have 
been for three dimensions. The Apollonian packing has 
dp = 2.4739 [68], and the values for space-filling bearings 
are again larger [59]. 

The fractal dimension is of particular interest here, be- 
cause it provides a means to characterize the properties 
of the disk areas. Melzak introduced the exponent of 
a packing, d r , defining it as the minimum value of s for 
which J2^Li r i 110 longer diverges. It was proved by Boyd 
that for the Apollonian packing that dp = d r [65]. To 
examine the divergence properties of this sum we can re- 
place the sum by an integral, because the divergence is 
controlled by the disks with small radii, for which the dis- 
tribution of radii is quasicontinuous. As this distribution 



follows a power law, p(r) 



[69], we have 



r s p(r)dr <■ 
(i ./o 

s+l— "I r ™ ox 



s + l-l3\ 
oo if s < - 1 



(42) 



Hence, dp = d r = (3 — 1. This allows the fractal dimen- 
sion to be determined from the exponent of the numer- 
ically obtained p(r) [66]. It follows that the area distri- 
bution is given by p(A) = A~( 1+dF/>2 * > for the 2D case 
and more generally for packings of d-dimensional hyper- 
spheres the volume distribution 



p(V) = v- {1+dF/d) « V~ 2 for large d. 



(43) 



Given the bounds for dp the exponent must lie between 
2 - 1/d and 2. 

One of the areas that is of most interest to us, and 
which is particularly relevant to the energy landscape 
networks, is the connection between the spatial proper- 
ties of the Apollonian packings and the topological prop- 
erties of the Apollonian networks. In Figure 10 we show 
the correlation between the disk area and degree. As ex- 
pected, the larger disks generally have a larger degree. 
However, for a given k there is a wide variety of disk ar- 
eas. The largest disks are associated with the crevices 
between the initial disks, whereas the smallest disks are 
obtained by following a spiral pathway in the network 
where each disk along the path is connected to one circle 
in each of the three previous generations. 

More specifically, the logarithmic average of the disk 
area for a given k closely follows a power-law. Assuming 
A(k) ~ k a and using the identity p(A)dA/dk = p(k) one 
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FIG. 10: The correlation of disk area with degree for the 
2D Apollonian packing (iterated to 15 generations). There is 
a dot representing each disk. The two lines with data points 
represent the average value of the area for a given degree, and 
correspond to both the normal and logarithmic averages. In 
addition lines corresponding to (fc/fc max ) 2 ' 62 and (fc/fc max ) 2 ' 428 
have been plotted for comparison. 
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FIG. 11: (colour online) The average area of the disks that 
a touch a disk against the area of that disk. The lines are 
binned averages for the Apollonian networks with t = 8 and 
a randomized version of that network. There is a data point 
for each disk in the Apollonian packing. Those represented 
by a cross are in contact with at least one of the three initial 
disks, and those represented by a diamond are not. 



can show that a = 2(7 - l)/(/3 - 1). For the 2D Apol- 
lonian network this leads to the prediction a = 2.428. A 
line with this exponent is plotted for comparison in Fig. 
10, and broadly follows the average (A(k)). By contrast, 
the average exp((logA(fc)}) has an exponent of 2.62. 

Although the correlations associated with the degree 
are most commonly studied, one can use Eq. (26) to de- 
fine an assortativity coefficient with respect to any prop- 
erty. Here, we examine the correlations in the areas 
of touching disks. From the behaviour of ta in Fig. 5 
one can see that there is some slight disassortativity, but 
much weaker than for the degree and with little differ- 
ence from that for the randomized graphs. However, if 
there was no variation in the areas of the disks with a 
given fc, (i.e. all the areas exactly obeyed the power-law 
dependence of k that characterizes the average) then ta 
and rfe would be very similar. Instead, because there 
is such a large scatter around the average values — over 
seven decades for the k = 3 disks in Fig. 10 — the effect 
of the lack of connections between nodes with the same 
degree only weakly carries over to disks with similar ar- 
eas. 

These effects can be examined in more detail by cal- 
culating A nn , the average area of the neighbours of a 
disk. Weak disassortativity is evident over the majority 
of the range of areas, except for the smallest disks which 
show strong assortativity (Figure 11). This latter effect 
is simply because the smaller disks in the last generation 
are connected to the smaller disks in the previous gen- 
erations, as with the spiral pathways mentioned above. 
Interestingly, Fig. 11 clearly divides the disks into two 
sets depending on whether they are in contact with one 
of the initial disks, and this is the source of the large 



fluctuations in the average value of A nn at intermediate 
values of the disk area. 



IV. DISCUSSION 

Our main motivation for studying the properties of the 
Apollonian networks is their potential to act as a useful 
model for the energy landscape networks. However, there 
is one major difference between the two systems. The 
Apollonian packings contain an infinite number of disks 
or hyperspheres, whereas configuration space is divided 
up into a finite number of basins, albeit a number that 
is an exponentially increasing function of the number of 
atoms in the system [70, 71]. 

There are two ways of creating a finite network from 
the complete Apollonian network. The first is to con- 
sider the network produced after a finite number of gen- 
erations, and is the one we have mainly used so far. The 
second is to consider the network containing only disks 
that are larger than a certain size. The wide distribution 
of areas for a given k in Fig. 10 indicates that their could 
potentially be significant differences. We know that the 
first will have a scale-free degree distribution, and the 
second a power-law distribution of radii up to their re- 
spective cutoffs, but what about the other way round. 

In Figure 12(a) the distribution of radii is shown af- 
ter different numbers of generations. These distributions 
approximately follow the expected power law for inter- 
mediate values of the radii, but this range becomes in- 
creasingly small as t decreases. Furthermore, at small r 
the lines curve away from this power law, because the 
finite packings only contain a small fraction of the total 
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FIG. 12: (a) The cumulative size distribution of the disks 
for an Apollonian packing after a finite number of genera- 
tions. The lines are labelled by their value of t, and an addi- 
tional straight line with the exponent expected for the com- 
plete Apollonian packing has been added for comparison, (b) 
The cumulative degree distribution for a series of Apollonian 
network where only disks with radii at least x% of that of 
the largest disk contribute. The lines are labelled by their 
value of x, and additional straight lines have been added for 
comparison, one of which has the exponent expected for the 
complete Apollonian packing. 



the small clusters used to generate the energy landscape 
networks [72] look quite like Fig. 12(a) suggesting that 
Apollonian networks with a given number of genera- 
tions are the more appropriate finite version for com- 
parison with these systems. Furthermore, there are then 
some useful parallels between t and N, the number of 
atoms in the cluster. For example, the number of min- 
ima increases exponentially with N and the number of 
disks/hypcrspheres in the Apollonian networks have a 
similar dependence on t (Eqs. 2 and 30). Similarly, as 
either t or JV increase, both types of networks become 
increasingly sparse (Eq. 5), have a smaller absolute value 
for the clustering coefficient C\ (Eq. 17), but a larger 
value relative to that for an Erdos-Renyi random graph 
(Eq. 20) [28, 29]. 

Other similarities between the two types of network 
include features that are quite common for scale-free 
networks, such as the dependence of vertex betweenness 
and local clustering coefficient on degree. Both are also 
disassortative [29], however, there is greater community 
structure in the Apollonian networks [72] . There are also 
similar relationships between the topological and spatial 
properties, such as for the dependence of disk or basin 
areas on k [72]. 

One of the interesting possibilities raised by the cur- 
rent study is the signature of the scale-free topology of 
the Apollonian network in the power-law behaviour of the 
disk areas. Currently, mapping out the whole network of 
connections between minima on an energy landscape is 
only feasible for systems of very small size. Neither are 
there methods available to construct a statistical rep- 
resentation of the whole network from a finite sample. 
Therefore, it is hard to test how generic is the scale-free 
behaviour observed for the clusters. However, the distri- 
bution of the hyperareas of the basins of attraction on 
a energy landscape is a static quantity that could po- 
tentially be statistically sampled for a system of arbi- 
trary size [32]. If this distribution exhibited a power law 
with exponent —2, (Eq. (43)), it would strongly suggest 
that underlying this was a scale-free energy landscape 
network, as for the Lennard- Jones clusters. Preliminary 
calculations indicate that this is the case [72]. 



number of disks in the complete packing with that r. 

The degree distributions for networks generated using 
a size cutoff are shown in Fig. 12(b). The distributions 
still follow a power-law, and are actually smoother, since 
k is no longer just restricted to the values given by Eq. 8. 
However, the exponent is slightly smaller than predicted 
by Eq. (13). The effect of the size cutoff is to only include 
the larger disks from the later generations, which are in 
turn more likely to be connected to the larger higher 
degree disks. For example, for a radius cutoff at 0.0001% 
of that of the largest disk, the first disks below the size 
cutoff occurred in the 13 th generation, and the last disks 
included were in the 706 th generation. 

Preliminary results for the basin area distributions for 



V. CONCLUSIONS 

We have analysed the properties of the contact net- 
works of space- filling packings of disks and hyperspheres, 
focussing on the Apollonian packing of two-dimensional 
disks. Their topological properties include a scale-free 
degree distribution whose exponent depends on the na- 
ture of the packing, high overall clustering, a local clus- 
tering coefficient that is inversely proportional to degree, 
disassortativity by degree and strong community struc- 
ture. 

These networks have many similarities to other deter- 
ministic scale-free networks introduced and analysed re- 
cently [36-38, 73, 74], but with the additional feature 
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that they have a well-defined spatial embedding. For this 
reason, we have suggested that these packings provide a 
useful model spatial scale-free network that may help to 
explain the properties of energy landscapes and the asso- 
ciated scale-free network of connected minima. In partic- 
ular, the scale-free topology of the Apollonian networks 



reflects the power-law distribution of disk sizes. Simi- 
larly, configuration space can be divided up into basins 
of attraction that surround the minima on the energy 
landscape. A similar power-law distribution for the hy- 
perareas of these basins might thus provide an explana- 
tion for the pattern of connections between the minima. 
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